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Einstein introduced the cosmological constant (A) as a fundamental parameter residing on the 
Left Hand Side (LHS) of his field equations for obtaining a finite density (p) for his model of static 
universe ( Einstein's static universe :ESU). For the sake of physical insight, we derive here the metric 
for ESU directly from Einstein equation, i.e., by considering both Einstein tensor Gik and energy 
momentum tensor T^. We find that in order that the fluid pressure and acceleration are uniform 
and finite despite the presence of a coordinate singularity, the effective density p e — p + A/8-7T = 0, 
where A is the cosmological constant. Under weak energy condition, this would imply p — A — for 
ESU. This means that if one would need to invoke a source of "repulsive gravity" in some model, 
(i) the model must be non-static, (ii) the repulsive gravity must be due to a "quintessence" or a 
"dark energy" fluid with negative pressure and appear on the right hand side (RHS) of the Einstein 
equation through Tij rather than through a fundamental constant (A) residing on the LHS of the 
same equation, and (iii) energy density of both normal matter and the "dark energy fluid" should 
be time dependent. In fact, the repulsive gravity would be due to a time independent fundamental 
constant A, it would be extremely difficult to understand why the associated energy density should 
be approximately 10 120 times lower than the value predicted by quantum gravity. On the other 
hand, for a dark energy fluid whose energy density is time dependent, it would be much easier to 
understand such an extremely low present energy density: the original initial value of the energy 
density of the fluid could be equal to the quantum gravity value while the present low value is due 
to decay with time. 
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INTROUDUCTION 

The "cosmological principle" demands that center of 
the universe lies everywhere and any fundamental ob- 
server must see the universe as isotropic. This means 
that, the space must be "maximally symmetric" [1] which 
demands that at each moment the metric of the space is 
same at all points and all directions [2]. By such con- 
siderations one can directly obtain not only the met- 
ric for Einstein's static universe (ESU) but the non- 
static Friedmann- Robertson- Walker (FRW) metric as 
well[l, 2] without using Einstein's equation 

Gik = 87T T ik (1) 

at all (G = c = 1). The Einstein tensor Gik, on the 
left hand side (LHS) of this equation, comprises only ge- 
ometric quantities, the structure of space time gik- On 
the other hand, the right hand side (RHS) of the same 
equation contains much more tangible quantity, Tj/., the 
energy momentum tensor of matter generating the space 
time structure gik or Gik- Thus any constraint imposed 
on gik from symmetry or other considerations must take 
its toll on the admissible form of as well. This means 
that constraints imposed on the LHS of Einstein equation 
must be reflected on the equation of state (EOS) of the 
fluid generating Gik- This latter physics aspect however 
may remain hidden in a purely symmetry driven mathe- 
matical derivation of cosmological metric which does not 
invoke Einstein equation at all. And in this paper, we 
would like to explore the tacit constraints imposed on 



by the requirements of the cosmological principle, "there 
is no preferred center", or "center lies everywhere". All 
the results obtained during this process will be also rati- 
fied by physical considerations as well. 

In effect we shall be dealing with spherically symmetric 
static solutions of Einstein equation with the inclusion of 
a cosmological constant A, a problem which was prob- 
ably first dealt by Weyl[3]. Since then several authors 
have considered such solutions and one particular case 
is due to Tolman[4, 5]. In recent times several authors 
have shown keen interest in this important problem [6, 7]. 
But, most comprehensive and elaborate treatment of this 
topic was due to Bohmer[8]. Since then Bohmer and his 
coworkers have tried to find even newer aspects of this 
problem which was earlier thought to be almost a dead 
topic [9, 10]. However such detail works mostly focused 
attention on isolated objects or stars where there (i) must 
be a pressure gradient, (ii) in general, a density gradient, 
(iii) a natural boundary R = Rb where fluid pressure 
must vanish, p(Rb) = and (iv) where a discuntinuity 
in the fluid density p might also occur. Further in such 
generalized studies, A was not considered as a fixed fun- 
damental constant contrary to what was in Einstein's 
mind when it was first introduced. 

In contrast, the present study, A = h.Einstein is a ba- 
sic constant and there is no question of modifying the 
value of A according to the needs of a general fluid so- 
lution. And of course, there must not be any natural 
"boundary" at all for the universe so that Copernican 
principle remains valid. It will be found that both the 
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metric for a star and ESU have a "coordinate singularity 
at R = a, where R is the luminosity distance and a is an 
appropriate integration constant. This singularity gets 
propagated in the hydrostatic balance equation as well. 
And in order to ensure that isotropic pressure gradient p' 
does not blow up at this singularity, relativistic stars are 
constructed so that Rb > a. Accordingly one can almost 
forget this singularity while studying static stars. But for 
the ESU one has to live with this singularity. And in this 
paper, we would study the behavior of physical quanti- 
ties like pressure and acceleration despite the presence of 
this singularity. 



FORMULATION 

We start with the assumption of spherical symmetry 
and consider a general form a static metric [2] 

ds 2 = e v dt 2 - e x dR 2 - i?W (2) 

where dfl 2 = d6 2 +sin 2 9d(j) 2 and R is the area coordinate. 
In particular, in a static universe, the luminosity distance 
turns out to be exactly R[3] 

d L = R (3) 

We also assume the cosmic fluid to be perfect with 

T ik = {p + p)uiU k - pg ik (4) 

where, p is the fluid density, p is the isotropic pressure, 
and u l is fluid 4- velocity. As is clear from Eq.(l), at the 
beginning, we do not consider any cosmological constant 
A; and the q component of the Einstein equation reads 

^ = ^=- e - x ih- X R) + h < 5 > 



This can be integrated to yield 



e- A = l- 



a(R) 
R 



where 



r 1 

a(R) = / 
Jo 



8irpR 2 dR 



(6) 



(7) 



Here the condition a(0) = has been used to ensure that 
e A is regular at R = 0. Thus whether, it is the interior 
solution of a star or the static universe, general form of 
the metric is 



ds 2 = e v dt 2 



dR 2 



l-a/R 



R 2 dQ 2 



(8) 



Note, as of now, we need not necessarily interpret a(R) as 
related to observed mass though what we have done looks 
like finding the interior of a relativistic star by working 



out the Schwarzschild interior solution. The important 
difference, however, is that while for a star there is a 
unique center and a boundary where the density may be 
discontinuous, for the universe, center is everywhere and 
there is no boundary, no exterior solution, and no den- 
sity discontinuity. Further even for a constant density 
star, there must be a pressure gradient because there is 
a unique center. But for the universe pressure gradient 
must vanish everywhere because the the center lies every- 
where and there is no boundary discontinuity. In fact for 
an isotropic and homogeneous universe, not only p and p, 
but all physically meaningful quantities, all geometrical 
scalars must the same everywhere. Such requirements 
may demand severe restriction on the admissible equa- 
tion of state (EOS) of the cosmic fluid. In contrast for a 
star, at best only p can be uniform in which case there 
would be a discontinuity at its boundary R = Rjj. 

Now we introduce the condition that the fluid must 
be homogeneous: p = uniform and thus both for the 
universe as well as a constant density star, one has 



a = — pR i 
3 1 

and the metric assumes the form 

dR 2 



ds z 



e v dt 2 - 



1 - {8n/3)pR 2 
If we introduce a parameter 



S = 



8np 



(9) 



R 2 dn 2 (10) 



(11) 



Eq.(10) would acquire the form 

dR 2 



ds 2 = e v dt 2 



1 - R 2 /S 2 

Further, if we introduce a new coordinate 
r = R/S 

we can rewrite the above metric as 



R 2 dn 2 (12) 



ds 2 = e v dt 2 - S 2 



dr 2 



1 - 7" 



+ r 2 d& 



(13) 



(14) 



Hence, the spatial section of both the universe and the 
interior of a constant density star is that of a 3-sphcrc, 
a space of constant curvature, a fact noted by Weyl[4]. 
This becomes clearer if we express 

r = R/S = sinx (15) 

to rewrite Eq.(14) as 

ds 2 = e v dt 2 - S 2 (d X 2 + sin 2 X dfl 2 ) (16) 

The coordinate singularity at R = S or a/R = 1 or r = 1 
now lies at the equator \ — 7r /2- This indicates that the 
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3-space of the problem occupies one-half of a 3-sphere of 
radius S. However such mathematical manipulation does 
not introduce any change of physics in the problem and 
one needs to be concerned with the singularity at R = a, 
or R = S, or, r = 1 or \ = 7r /2- 

Let us suppose that somehow the density of the star 
would be reduced to zero everywhere. In such a case, the 
space time must become flat with a metric 



and an effective pressure 



ds 2 = dt 2 - dR 2 - R 2 dil 2 



(17) 



From the viewpoint of metric(8) it would immediately be 
apparent that one should then have e v = 1 and a/R = 0. 
But from the viewpoint of metric(14), it would appear 
that in such a case one would have S — oo, r — and 
R = Sr > 0. Why would we have r = in the latter 
case? This is so because the general form of a space time 
with constant spatial curvature is 



ds 2 



° v dt 2 



dr 2 



1-Kr 2 



dn- 



(18) 



where K = +1, 0, —1 corresponds to closed, flat and open 
space time respectively. And once we presume K = +1 
and adopt metric(14), we presume p > and then it 
becomes difficult to have a smooth transition to a p — > 
case. Hence, if one would try to arrive at a flat space 
(K = 0) by imagining that it is due to a closed space 
of infinite radius, one would land up with the condition 
r = 0. In fact, if in Eq.(4), had we taken, 



T ° = Kp 



(19) 



instead of p, we would indeed have obtained Eq.(18) di- 
rectly in lieu of Eq.(14)- Then we would have found that 
the condition for flatness is K = 0. But when we do 
take Tq = p and implicitly assume p > 0, we presume 
K = 1 and exclude the possibility that one might have 
K = too. In such a case, the condition for flatness 
would appear as r = as mentioned above. 



Cosmological Constant 

It is well known that in order to obtain a static uni- 
verse which would be closed and finite radius, Einstein 
modified Eq.(l) into 



G ik + Ag ik = 8w T ik 



(20) 



apparently implying that either A is a fundamental con- 
stant like 87rG/c 4 or a basic scalar, like the Ricci scalar 
1Z appearing within Gik- It follows that, now one has an 
"effective" density [8] 



P + 



A 
8^ 



(21) 



A 



So instead of Eq.(7), one now obtains 

_ A _ a e (R) 
R 



where 



fit 



a e = 8np e R 2 dR 
Jo 

One also finds 

ds 2 = e v dt 2 



dR 2 



R 2 dil 2 



1 - Ote/R 

And for a constant density case, one has 

dR 2 



(22) 



(23) 



(24) 



(25) 



* 2 = e ^- l-(8,/3WP + * dn2 (26) 



and 



8irp e 



(27) 



Eqs. (12-16) remains unaltered in the presence of A. Thus 
now, if the space time has to be flat, in addition toe" = 1, 
one must have 



(28) 



This shows that if the original definition of "vacuum" is 
p = 0, as if, in the presence of a A, it gets modified to 
Pe = 0. 



ACCELERATION SCALAR AND SINGULARITY 

The fact that for a spherically symmetric static sys- 
tem, one should have a/R < 1 has already been 
investigated[ll, 12]. The basic reason for this is not dif- 
ficult to see. In general, static or non-static, for spher- 
ically symmetric space time occurrence of a(r,t)/R > 
corresponds to a formation of a "trapped surface" and 
the condition a(r,t)/R = 1 marks the formation of an 
"apparent horizon". Though for a non-static system, it 
might be possible to have trapped surfaces or horizons, 
for a static system they are not allowed. This is so be- 
cause once a trapped surface would be there, stellar mat- 
ter would be inexorably pulled towards the central point 
of symmetry and thus matter would soon end up in a 
point singularity rather than as an extended static ob- 
ject. While this is definitely not allowed for a static star, 
this problem would be much more severe for cosmology 
because "center of symmetry lies everywhere" . 
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Following the case of a Schwarzschild black hole space 
time, generally, it is believed that this R = a singularity 
is a mere coordinate one even in the presence of matter. 
But as we would see below, the coordinate independent 
scalar acceleration of the fluid would blow up unless se- 
vere constraints are imposed on the fluid EOS. 

It may be recalled that Einstein was very much con- 
cerned about this R = a singularity and constructed 
a static model of a fluid where test particles arc mov- 
ing in randomly oriented circular orbits under their own 
gravitational field. While the radial stress of the fluid is 
zero, tangential stresses are finite. This configuaration is 
known as "Einstein Cluster" and Einstein showed that, 
the speed of the orbiting particles would be equal to the 
speed of light c at R = a [13]. Thus, he pleaded that 
there cannot be any R = a or R < a situation. Later 
it turned out that Einstein's cluster indeed corresponded 
to some well defined interior Schwarzschild solution stud- 
ied by Florides[14, 15]. And now Bohmcr & Lobo[16] 
have shown that Einstein's intuition was correct, and the 
R = a singularity is indeed a curvature singularity. This 
however does not at all mean that all R = a singulari- 
ties are curvature singularities [8]. In fact, in the present 
static case of isotropic pressure, it would be found that 
they should be regions with zero curvature singularity to 
avoid an acceleration or pressure gradient singularities. 

For any static spherically symmetric fluid, one can eas- 
ily find the acceleration[10] 



(29) 



In spherical symmetry, only one component of accelera- 
tion survives 



er x v' 



(30) 



Using Eq.(33), we obtain 



a + 8ttR 3 p 
2RS/Y = c7J% 



(35) 



If A would be included, then the acceleration scalar would 
be given by 



a e + 8irR 3 p e 



2R 2 y / l-a e /R 



(36) 



It is now clear that if a static fluid would tread upon 
the R = a e singularity, its acceleration could be infinite 
and it cannot be at rest. And this is the physical reason 
that for a static fluid one must have R > a e everywhere. 
And in case, the manifold would cover R — a e , one must 
satisfy 



a e + 8itR p e = 0; R = a e 



(37) 



in an attempt to keep a regular. Using R — a e in the 
foregoing equation, we have 

a e (l + 8na 2 ePe ) = (38) 

i.e., one must have cither 

a e = (39) 



or, 



Pe{R = ae) = 



(40) 



or both of the above two conditions. If we assume that 
the minimum value of p is zero and it cannot be negative, 
we will have p e = —A/8n. Then Eq.(40) would yield 



How to evaluate this a R l Again, for the sake of easy 
understanding, we first do not consider any A and write 
down the ^ component of Einstein equation [2], 



In view of Eq.(30), let us rewrite this equation as 
er x a + 8irpR 3 



R 



R 3 



or 



o a + 8ttR 3 p 
a = 



(31) 



(32) 



(33) 



Since there is only one component of a 1 , the scalar accel- 
eration becomes 



a = yf-tfai = a R yJ\g RB \ 



(34) 



R = S 



(41) 



But one is still not sure whether additional condition such 
as Eq.(39) is needed to really ensure that a is indeed finite 
at R — a e . For a fluid with p e > 0, the safest way to 
avoid the R — a e singularity will simply be to ensure that 
R> a e . The fluid can do so by choosing an appropriate 
density profile and by ensuring that its outer boundary 

Rb > a e (42) 

For a constant density star or the ESU, we have 

R(p e + 3p e ) 



a = k- 



yjl - 2np e R 2 
where k = 4tt/3. In terms of r, we obtain 

,r(p e + 3p e ) 



a = kS- 



(43) 



(44) 



■5 



It is clear that the sufficient condition for avoiding the 
R = a e singularity in this case is 



R b > a; R< S 

In terms of density, this means 

3 

Pe > 



8nRl 



(45) 



(46) 



Thus, if one would imagine a region with p = 0, i.e., 
p e = A/8-7T, one should restrict the radius of this vacuum 

as 



R b < 



(47) 



And in case, one would have p — A/Att, then the above 
restriction would become 



R b < 



(48) 



and which cannot be satisfied [8]. Thus for a constant 
density star with R^ = y/T/X, one must have p > A/4-7T. 
And in case this condition would be violated, one must 
critically analyze the additional constraint on the EOS 
which would prevent p' e from blowing up at R — a e . It 
may be seen that the p = A/ An EOS corresponds to 



if p = 0. 



p e + 3p e = 



TOV EQUATION 



(49) 



Recall that local energy momentum conservation equa 
in T*. k = immediately leads to [2] 

, -v 



p + p 



(50) 



Further if we again consider the ^ component of Einstein 
Eq. and combine it with Eq.(50), we will obtain the TOV 
equation for hydrostatic balance for any self-gravitating 
static fluid: 



, (p + p)(a + 8TrR 3 p) 
P =-- 



2R 2 {1 -a/R) 



(51) 



Since the effect of A gets included if one replaces p and 
p by their "effective" values [8] 



, = _ (Pe + Pe)(a e + &TTR 3 Pe) 



2R 2 (l~a e /R) 



(52) 



Further since p' — p' el in a complete symmetric manner, 
we rewrite the above equation as 



Pe = 



(Pe +Pe){a e + 8irR 3 p e ) 
2R 2 (l-a e /R) 



(53) 



This equation strongly suggests that if the definition of a 
"dust" in the absence of A is p = 0, as if, in the presence 
of A, dust EOS should be p e — 0. Similarly, if the original 
EOS of "vacuum" is p = 0, in the presence of A, the 
vacuum EOS is p e = 0, a hint we have already found. 

Note, even now, there is no need to interpret a in terms 
of any exterior boundary condition, and thus TOV equa- 
tion is valid in any spherically symmetric static GR prob- 
lem including ESU. By using Eq.(36), it is interesting to 
rewrite the TOV Eq. in terms of the acceleration scalar 



Pe 



a{Pe + Pe) 

y/1 -a e /R 



(54) 



Clearly, there is a singularity in the denominator of TOV 
Eq. at R = a e . If we write 



x 



= y/1 - a e /R 



(55) 



it is seen that while the singularity in acceleration a <~ 
x , for the pressure gradient it is much stronger p' ~ 
x~ 2 . And if the fluid would cover the R = a e singularity, 
regularity of a may not be sufficient to ensure regularity 
of p' . In fact Eq.(54) suggests that one might require the 
additional constraint 



p e +p e = 0; R = S 



(56) 



to tame the p' singularity at R = S. 

Now let us use the condition p e — uniform which is 
equally valid for the interior solution of a constant density 
star as well as the ESU: 



p' = -kR 



(fie +Pe)(Pe + 3p e ) 



1 - 2np e R 2 

In terms of the normalized coordinate r, one finds 

j(p e +p e )(Pe + 3Pe) 



P ' = -kS- 



And in term of a, we have 



1-r 2 



, a(p e +p e ) 
P = ~ 



and 



(57) 



(58) 



(59) 



(60) 



y/1 - 2n Pe R 2 

, _ a{pe +p e ) 

p ~~ 

CONSTANT DENSITY STAR 



For the interior solution of a constant density relativis- 
tic star one, one finds[8] 



P c + P 
P + P 



(61) 
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where p c is the central pressure. This relation is inter- 
esting because it does not involve any exterior bounday 
condition. Also note that since 



or both 



P- 



we can rewrite Eq.(61) as 



Pe + Pe 



Pe + Pe 
Pe + Pe 



(62) 



(63) 



Since for a star, there must be a pressure gradient, 
p c 7^ p(R), and one really cannot reduce e v (R) to a con- 
stant value independent of R. But if, from the mathemat- 
ical viewpoint, one would still demand that it should be 
possible to set up a time orthogonal Gaussian coordinate 
system where dt 2 = dr 2 , one would unknowingly kill the 
pressure gradient and set p c = p in Eq.(61). In the ab- 
sence of A, for a static configuration, p — if p = 0, and 
the star would vanish under the assumption of e v = 1! 

When A is present, the expression for pressure for a 
constant density star is [8] 



p(R) 



A/47T - 1 + C cos x 
3 — cos x 



where 



C = 



3p C +p- A/47T _ Zp C e + Pe 



Pe + Pe 



(64) 



(65) 



Now it transpires that even when p = to honor 
e v = 1, there would still be a finite density p = K/A-k. 
But one attains this finite density at a huge cost because 
it now turns out that the boundary of the star merges 
with the coordinate singularity R b — S where the metric 
behaves badly. For a real star, this singularity can al- 
ways be avoided by allowing e v > 1 which in turn means 
by having p > A/ Air. Thus a "world time cannot be 
introduced for a general relativistic star having a finite 
density. 



Verification From TOV Equation 

Although the case of constant density stars in the pres- 
ence of cosmological constants has been studied in great 
detail[8], to our knowledge, nobody actually verified the 
validity of the solutions in case one would tread on the 
coordinate singularity at R — S. Since the denominator 
of Eq.(53) becomes zero at R = S, in order that p' docs 
not blow up there, one must have either 



Pe+Pe = 0; 



R = S 



(66) 



Pe + Pe = Pe + 'iPe = 0; R= S 



(68) 



If the last constraint would be satisfied, one would imme- 
diately obtain p e = p e = 0. However, it is possible that 
Eq.(68) is not satisfied and only one of the less rigourous 
conditions(66) or (67) is satisfied. In particular, to avoid 
occurrence of negative pressure, let us assume that only 
Eq.(67) is satisfied at R — S. Since p' attains a 0/0 form 
at R = S, let us study the nature of p' at this singularity 
by using L' Hospital rule. For this, let us first write 



Pe = 



f(R) 
g{R) 



where 



and 



f(R) = K R(p e + Pe)(Pe + 3p e ) 



g(R) = 2 K p e R 2 - 1 
Since p' = 0, we find that 



(69) 

(70) 
(71) 



/' = n(p e +p e )(p e + 3p e ) + nRp e [3(p e +p e ) + (p e + 3p e )] 

(72) 

Since, we have already considered p e + 3p e = 0, we obtain 
a reduced expression for 



f = 3 K Rp' e (p e +p e ) 
On the other hand, 

g' = 4np e R 

So that 

4=3^ 

9 4p e 



(73) 
(74) 

(75) 



And at, the singularity, R = S 7 by 1' Hospital rule, we 
obtain 



p' e = lim 

R^S g' 



(76) 



From Eqs.(75) and (76) we obtain the required condi- 
tion: 



i.e., 



HPe +Pe) = 4p e 



p e (R = s) = 3 Pe (R = S) 



(77) 



(78) 



or 



which looks like the EOS of incoherent radiation! By 
combining Eqs.(67) and (78), it becomes clear that in 
order that p' does not blow up at the "coordinate singu- 
larity" at R — S, the constant density star must have 



p e + 'ip e — 0; R = S 



(67) 



Pe(R=S)= P e(R=S)=0 



(79) 
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Thus instead of p = we would obtain p e = at 
R = S. This means that if indeed A > 0, the solu- 
tion must avoid R = S singularity. More importantly, 
since p e = constant, we find that, if the solution would 
indeed extend upto R = 5, we must have p e = 0. Thus 
all finite density stars must avoid the a e /R= 1 orr = f 
singularity. 



STATIC UNIVERSE 

In cosmological case, there must be a universal time 
which would be the proper time of all fundamental ob- 
servers, i.e., dt 2 = dr 2 . This demands that one must be 
able to set e v — 1 so that, the Eq.(14) would become the 
ESU metric 



ds 2 = dt 2 



dr 2 



l-r 2 



dfr 



(80) 



Note that the value of S for ESU is still given by 
Eq.(27)[10] and thus one is justified in deducing ESU 
metric indeed by setting e v = 1. Since the Einstein equa- 
tion tells that the metric is essentially determined by , 
such an important change of setting e" = 1 on the LHS 
of this equation must be endorsed by the RHS, i.e., by 
the admissible forms of fluid EOS. What are those condi- 
tions? To explore them, we note from Eqs.(50) and (51) 
that the most general form of static, spherically symmet- 
ric Einstein Eq. yields 



, a e + 8irR 3 p e 
~ R 2 {1 -a e /R) 

In constant density case, this condition becomes 



v' = -2k 



R(p e + 3p e ) 



1 - 2np e R 2 
And in term of r, this condition becomes 



(81) 



(82) 



(83) 



In order that c v — 1, one must have v' — 0, and from the 
foregoing equation, it is immediately clear that then we 
must have atleast 



p e + 3p e = 



(84) 



In any case, we must have p' e = everywhere including 
R = S or r = 1. Then if we follow the L' Hospital 
treatment of the previous section, we would find, we must 
have 



Pe = 



(85) 



everywhere because p e — constant. Although, we have 
already obtained this important result, it would be in- 
teresting to obtain this result from somewhat different 
routes. 



First consider Eq.(59) and write 

f(R) = a(p e +p e ) 

and 



so that 



and 



g(R) - -^\-2np e R 2 



f'=p' e a 



g' = 2np e (l - 2 K p e R 2 )- 1/2 
Now applying L'Hospital rule at R = S, we find 



, , a-v/l — 2np e R 2 
P = P ^ 5 

Using Eq.(43) in the above Eq., we find 

Pe + iPe 



2 Pe 



1 



implying 



Pe + Pe = 



(86) 
(87) 
(88) 
(89) 

(90) 

(91) 
(92) 



i.e., both the factors in the numerator of Eq.(58) should 
vanish to ensure regularity of p' at R = S. It may be 
also of some interest to extend this study by directly 
considering the radial variable as r rather than R. Note 
that 



, 1 dp e 



(93) 



Now if we denote differentiation by r with a prime, we 
rewrite Eq.(60) as 



, _ Sa(p e +p e ) 
P ~~ ' a/T - 72 



(94) 



so that 



and 



Then we have 



and 



f(r)=Sa(p e + Pe ) (95) 

g(r) = -\J\-r 2 (96) 

/' = Sap' (97) 

*' = 7CT (98) 
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Again applying L'Hospital rule for the limit of p' at R = 
S, we find 



SaVT^r 2 



= 1 



(99) 



Inserting the expression for a from Eq.(44), we obtain 
the interesting relation 



K S 2 (p e + 3p e ) = 1 



or 



3 



47r(/9 e + 3p e ) 



(100) 



(101) 



Comparison of Eqs.(27) and (101) would again convey 
one of the hidden mesages for the ESU as well as for 
a constant density star which is attemting to suppress 
its pressure gradient, namely 3p e = p e - The sum and 
substance of the entire exercise is that for the ESU, in 
oder that p' e is indeed zero everywhere, one must have 



■i Pe = 



(102) 



But unlike the constant density stellar case now we can- 
not escape confronting this singularity by demanding that 
p > A/47T. And hence we must accept the fact that for 
the ESU fluid p e — p e — everywhere. This means that 



p + A/47T = 0; p=-4A/47r 



(103) 



atlcast for the ESU. The weak energy condition demands 
that p > and thus we find that A = 0! 



EINSTEIN'S SOLUTION 

If one would directly use Eq.(80) into Einstein Eq.(20), 
one would be led to [5] 



4ir(p + 3p) = A 



and 



4ir(p + p) = 



S 2 



Einstein chose, p = EOS which led to 



(104) 



(105) 



(106) 



Note that Eqs.(104) and (105) may be rewritten as 

4ir( Pe + 3p e ) = (107) 

and 

1 



47r(p e +p e ) 



S 2 



(108) 



But, if we recall, Eq.(92), actually, p e + p e 
S = oo. Further Eqs.(100-101) showed that 



47r(/9 e + 3p e ) 



_3_ 

S 2 



so that 



(109) 



And a comparison of Eqs.(108) and (109) again shows 
that, S = 00 and p e = 0. As mentioned before, for weak 
energy condition, this would mean p = A = 0. Note that 
it is indeed possible to have a situation where p rest /e^> 
1, where e is the proper internal energy density. In such 
a case, one might approximately write p w 0. But this 
does not mean that pressure is strictly zero. A strict 
p = EOS is possible only when p = 0. In such a case, 
Eq.(104) would again yield A = 0. 

If we ignore such physical and mathematical regular- 
ity or self-consistency considerations, Einstein's solution 
would suggest A = A-nGp/c 2 . 

Further, if one would assume that (i) ESU correctly 
describes the real universe and (ii) the true mean den- 
sity of the real universe in its totality (about which we 
may never have absolute knowledge) is equal to the mean 
density of the patch of the observed universe, i.e., 



Ppatch — Ptrue P ' 



io~ 31 g cur 3 



one would obtain[3] 



A - 10 



" 58 cm 2 



(110) 



(111) 



Now ponder over the fact that until the advent of big op- 
tical telescopes in 1920 or so, for hundreds of years, most 
of the astronomers thought that total universe was noth- 
ing but our milkyway galaxy. And if would have been 
possible to measure the mean density of the Milkyway, 
one might have concluded that 



Ppatch — Ptrue P ' 



10~ 13 g cm" 3 



and accordingly 

A - 10~ 40 cm 2 before 1920 



(112) 



(113) 



As late as 1970, we hardly had any idea that the galaxy 
clusters and superclusters are actially distributed as "fila- 
ments" and "walls" around huge voids whose dimensions 
could be as large as ~ 280 Mpc[17]. The recent Sloan 
Digital Sky Survey has revealed structures of dimension 
<~ 500 Mpc[18]. Thus, we cannot rule out the possibility 
that that eventually, it might be found that the entire 
patch of presently observed universe also lies on the wall 
of a larger void. If the mean density of observed universe 
would be revised by future observations, would we again 
revise the value of the fundamental constant A? 

Many cosmologists believe that within the patch of the 
observed universe, galaxies are distributed in fractal like 
pattern with length scale I ~ 100 Mpc and a fractal di- 
mension D <~ 2 [19-21]. One may recall here a succint 
comment by Wu, Lahav & Rees[22] 
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"The universe is inhomogeneous - and essentially frac- 
tal - on the scale of galaxies and clusters of galaxies, but 
most cosmologists believe that on larger scales it becomes 
isotropic and homogeneous." 

In a fractal distribution, the number of galaxies in- 
creases with R as 

N(< R) ~ R D (114) 

and the number density of galaxies varies as 

n(< R) ~ R°- 3 (115) 

Thus, in principle, for a fractal distribution over suffi- 
ciently large scale, one may have n — > 0, p — > 0. Most 
of the fractal patterns however tend to have a preferred 
center and thus not suitable for an eventually homoge- 
neous and isotropic universe. But there could be some 
fractal structures like Levy Dust which may not have any 
preferred center. One may have a mean p = also for an 
appropriate infinite multi universe[23, 24]. Suppose the 
entire observed universe where we live is only a speckle 
in a multi universe comprising infinite such speckles. If 
these speckles or "cosmic atoms" are separated from one 
another by infinite voids, then the mean density of the 
infinitely diluted multi universe would be zero. 

At the cost of sounding repetitive, let us again raise the 
old question: if A would indeed be due to some quantum 
mechanical effect, why the value of A obtained under the 
assumption of p pa tch — Ptrue falls short of the theoretical 
value by an approximate factor of 10 120 ? And if the 
various field modes would cancel one another to generate 
a small A, why do they not cancel exactly to result A = 
0? In fact there are some theoretical estimates by which 
A = 0[25]. 

In physics, we have various fundamental constants like 
G, c, h, m e etc. But in no case there is any evidence 
that the values of such constants directly depend of am- 
bient factors like mean density of universe. In fact, all 
truly fundamental factors allow themselves to be deter- 
mined with high precision by judicious combination of 
theory and experiments. The only exception here is the 
supposed A! 

Thus if one would indeed need to invoke a source of 
"repulsive gravity" in some model, the source cannot be 
a fundamental time independent constant (A) sitting on 
the LHS of Einstein equation (20). On the other hand, 
such repulsive gravity could be due to a "dark energy" 
fluid represented by an appropriate sitting on the 
RHS of Einstein equation. In such a case, in a non- 
static case, the energy density of the dark energy fluid 
also must be non-static and except for a fiducial initial 
value, such an energy density may not considered as any 
"fundamental parameter" . 



Faintness of Distant Supernovae 

Over the past decade, cosmologists have gathered ev- 
idences that the luminosities of distant Type la supcr- 
nove are lower than what is expected in the standard 
FRW model. Here, it is assumed that the Type la super- 
novae luminosities are indeed "standard candles" despite 
the likely evolution of luminosities due to likely cosmic 
evolution. Even if the Type la supernovae would be stan- 
dard candles without appreciable cosmic evolution, there 
is some chance that their faintness could be due to opaci- 
ties of distant inter galactic medium. For example, it has 
been pointed out that Lyman Alpha clouds might in- 
troduce some non-transparency for optical emission from 
distant universe [26]. It has been also argued that the at- 
mosphere of planets could be additional sources of opac- 
ities for very distant supernovae lights [27]. 

If we ignore this possibility of dimming of distant su- 
pernovae due to some intergalactic opacity, for the FRW 
model, one needs to postulate an "accelerating universe" 
where the acceleration is driven by some source of repul- 
sive gravity. But as discussed in the previous subsection, 
this repulsive gravity should be due to a dark energy fluid 
whose energy density is decayed from the quantum grav- 
ity value by a factor of almost 10 120 . And this repulsive 
source of gravity should appear from rather than as a 
fundamental constant from the LHS of Einstein equation. 

Also, one should leave here an open window for some 
hitherto unknown physics or surprise for understanding 
this new phenomenon. For instance even now, there is 
no firm evidence that Copernican principle, the basis of 
the standard cosmology is indeed valid[28]. Also, con- 
sider the fact that we are able to "see" galaxies and mea- 
sure their redshifts because they are radiating. But in 
the standard cosmologies, galaxies are considered as non- 
radiating neutral "dust" particles. 

SUMMARY AND DISCUSSIONS 

Both a spherically symmetric static star and the ESU 
result from the same spherically symmetrical form of Ein- 
stein equations. If j| component of the Einstein equation 
would be studied, one would obtain, the expression for 
acceleration scalar and condition for hydrostatic balance 
in both the cases. In fact, Bohmer[8] did obtain one of the 
constraints to be followed by the ESU fluid by demanding 
that the pressure gradient p' must vanish for the ESU: 
p e + 'ip e = 0. We found that the singularity appearing in 
the (R, R) component of the metric tensor in all spher- 
ically symmetric static configurations do propage in the 
expressions for acceleration scalar a and pressure gradi- 
ent p' . While, for the former, the nature of singularity is 
a <~ a; -1 , for the latter it is much stronger p' ~ x~ 2 . In 
fact, the requirement p e + 3p e = is the basic condition 
for ensuring that the a does not blow up at x = or r = 1 . 
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This basic condition however need not be the sufficient 
condition for ensuring that a is finite at x — 0. Further, 
since p' tends to blow up much faster at x = 0, one must 
require an additional condition to ensure regularity of p'. 
We found that the latter requirement is p e = 3p e . This 
means that if the constant density static star would have 
Rb = S, one must have p e = 0. But one can always 
avoid this singularity by ensuring that p > A/4w. How- 
ever, previous detail studies of perfect fluid spheres were 
primarily intended for isolated systems and not for the 
universe, and to our knowledge, nobody tried to study 
the remedy for either the ~ x" 1 or ~ x~ 2 singularities 
at x = 0. 

Our result that whether it is the ESU or a star with 
Rb = S should have p e = to tame the p' <~ x~ 2 singu- 
larity can be obtained probably from other directions as 
well. Consider the interior metric(7) where t is the time 
measured by a clock at a fixed R. In the exterior vacuum 
Schwarzschild metric 



ds 2 = (1 - a e /R)dt 2 



dR 2 



1 — a e /R 



- R 2 dQ 2 (116) 



too, time is measured by a clock resting at R = R. 
Suppose we match the interior and exterior solutions at 
R = Rb- If such a matching is possible, we would have 



e v » = 1 - a e /R 



(117) 



In particular, consider the constant density solution[8] 
where 



p c + p 



(118) 



If one would have a = (8/9) Rb, central pressure would 
blow up (assuming p > 0). For a uniform density star 
this borderline would correspond to 



8/9 



(122) 



But suppose there is a solution which appears to vio- 
late this constraint, i.e., r > 8/9. Then how to ensure 
that central pressure is still finite? The only solution out 
of this dilemma would be to set a = or p = 0. From 
this more stringent condition, it should be clear why Ein- 
stein's original static universe must have p = 0. 

Does the fate of the ESU improve after the incorpora- 
tion of a positive A? From Eq.(3.35) of ref.[8], it turns 
out that now the modified Buchadhal constraint is 



a. 2/ _ A 
R 9 V Airp 



(123) 



Thus with a positive A, this constraint gets even tighter. 
With p = A/47T, one obtains 



a e 2 
~R < 3 



(124) 



which means, for a constant density star, one must have 



r < 2/3 



(125) 



to ensure that pressure does not blow up at the center r = 
0! For the ESU, this would mean blowing up of pressure 
everywhere because "center would lie everywhere" . And 
the only solution to get rid of this problem is is accept 
the fact that a e = for ESU. 



In order that e Vb — > 1, it is clear that, one must have 

a e /R b = ^p e Rl - (119) 

which implies 

Pe -> (120) 



It is also seen from Eq.(118) that e Vb — > 1, only if p c — > 0. 
It so happens that, one fails to match the interior solution 
with the exterior vacuum solution(7) at p = A/ An. As 
we have found, one fails to do so because gnu blows up 
in Eqs.(10) and (116) in this case. 



Buchdahl Inequality 

It is well known that any spherically symmetric static 
configuration, homogeneous or inhomogeneous, satisfies 
a constraint much stronger than the Rb > a or a < 
Rb constraint. In the absence of a A, this is known as 
Buchdahl inequality [1, 26]: 



a 8 
Rb < 9 



(121) 



Mass Function 

Since for a fluid in an asympotictially flat space time, 
one can define a "mass function" in terms of which ac- 
celeration scalar 



M + AttR^p 
R 2 ^l - 2M/R 



(126) 



we can now identify a(R) in Eq.(36) as twice the quasilo- 
cal mass-energy of the fluid 



a{R) = 2M(R) 



(127) 



In the presence of a A, the quasilocal mass will be M e — 

Oe/2. 

In general relativity, the gravitational mass of a sta- 
tionary system is [4] 



M = /(T ° + i°) V=9d 3 x 



(128) 



where t\ is energy momentum pseudo tensor and g is the 
determinant of the metric tensor. One can work out t® 
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from the metric of the universe. Probably starting 
from Rosen [27], many authors have worked out the mass 
energy of a closed universe and all of them have concluded 
that for a closed universe M = [28-32]. This too would 
tell that, in the absence of A, p = 0, and, in general 
p e = 0. For the ESU, it is possible to confirm this result 
irrespective of the value of In all coordinate systems, 
for a static fluid, one obtains [28-30] 

M e = J (p+3p-A/4n) ^ d 3 x = J (p e +3p e ) J~g d 3 



(129) 



(130) 



In view of Eq.(84) we then directly obtain 
M e = 

Then from Eq.(24), we obtain p e = 
Poisson Equation 



It is known that the RHS of Poission equation indicates 
the source of gravity. For a spherically symmetric static 
fluid the Poission's equation is 

V 2 v / ffoo = 4:iry/g^(p + 3p- A/4tt) = An^/g^(p e + 3p e ) 



This shows that the source density of gravity is 

Pg = V500(pe + 3p e ) 

And when goo = 1, one will obtain 

Pg = (Pe + 3p e ) 



(131) 



(132) 



(133) 



In either case, the RHS of Eq.(131) is zero when p e + 
3p e = 0, which is the case for both a constant density 
star (if it would be extended up to R = S) or the ESU. 
In all such cases, one must necessarily have goo — c v — 
constant. And since source of gravity is zero, the space 
time must be flat. And a space time is flat when p e = 0, 
as obtained by us. Alternatively, if one would first focus 
attention on the LHS of Eq.(131), it seems that if one 
would set e v = constant, the source of gravity will vanish, 
again in which case, one should have p e = 0. 



A Simple Reason 

For a constant density case, we may rearrange the TOV 
equation(58) as 



-P 



kS 



(pe + Pe)(Pe + 3p e ) 1 - r 2 

We may also rewrite the acceleration equation(44) as 
a kS r 



(134) 



Pe + 3p e Vl - r 2 



(135) 



The RHS of Eqs.(134) and (135) obviously depends on 
r if r would indeed be a free parameter. For a constant 
density star, the LHS of the same equations too must 
depend on r. And they do depend on r because p e = 
p e (i~) even though p e = constant. But suppose we would 
like to freeze the r dependence of p to make e" = 1. 
Then the value of r on the RHS of Eqs.(134) and (135) 
too must be frozen. And since R = Sr > during the 
freezing process, the only solution here would be to adopt 



oo; 



0; 



R> 



(136) 



We have already discussed that occurrence of Eq.(136) 
actually signifies that the spatial section is flat: K = 0. 
This would have been more transparent had we taken 



n 



-Kp 



alongwith Eq.(19), i.e., had we written 

Tik = K[(p + p)uiU k + pg ik ] 



(137) 



(138) 



instead of Eq.(4). Further had we also taken the cosmo- 
logical constant as KA instead of A, Eqs.(134) and (135) 
would have appeared as 



and 



KkS r 



(pe+Pe)(Pe + 'iPe) 1 - Kr 2 



KkS 



Pe + 'iPe VI - Kr 2 



(139) 



(140) 



respectively. In such a case, it would have been seen 
immediately that the requirement that the RHS of 
Eqs.(139) and (140) are independent of r, one must have 
is K = 0. And since we have already set e" = 1, it 
means, the space time must be flat to ensure that the 
RHS of Eqs.(134) & (135) or (139) & (140) are indeed 
independent of r. Hence one must have Kp e = in such 
a case. So will be the case for ESU. 



CONCLUSIONS 

Cosmological constant was introduced as a fundamen- 
tal constant on the LHS of the Einstein equation. And 
we found that, if A would indeed be considered as a 
time independent fundamental constant, its value is ex- 
actly zero rather than lower than the quantum gravity 
value by a factor of ~ 10 120 . This, we believe, has been 
shown for the first time and is an important result in its 
own right. Thus the Einstein equation should be Eq.(l): 
Git — 87r To- rather than Eq.(20). However, in principle, 
one is free to choose any kind of source of gravity, i.e., 
Tjfe. In particular, one may choose T ik as a blend of a nor- 
mal fluid with positive pressure and a "dark energy fluid" 
with negative pressure. The latter fluid would then act as 
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a source of "repulsive gravity" . The most popular model 
to explaining the extra faintness of distant supernova in- 
vokes such quintessence or dark energy fluid. Since in a 
non-static universe model, all components of includ- 
ing the energy density of dark energy fluid would be time 
dependent, it would be easier to understand why the re- 
quired present day density of this dark energy fluid is 
~ 10 120 times lower than the quantum gravity value. On 
the other hand, had this "repulsive gravity" been due to 
a time independent A, this task would have been much 
more difficult. 
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